Pre-management of Non-local Decoherence 
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The possibility of pre-managing a non-local decoherence process with limited resources based on 
initial conditions is an open question. A generic two-qubit system that undergoes local but different 
amplitude-damping dissipations is studied to analyze the relation between entanglement sudden 
death (ESD) properties and initial conditions. ESD-free, ESD-tolerable and ESD-no-go phases are 
shown to be pre-manageable through the initial entanglement and excitation probabilities. Our 
results may provide useful guidance to state preparations in practice. 

PACS numbers: 03.65.Ud, 03.65.Yz, 42.50.Pq, 75.10.Jm 



BACKGROUND 

Coherence is one of the central requirements for the 
success of quantum computation and quantum commu- 
nication tasks [![. However, the existence of unavoid- 
able environmental interactions will induce both local 
and non-local decoherences In particular, this irre- 
versible process may lead to the finite-time disappearance 
of non-local coherence Q when local coherences still re- 
main non-zero. This is well known as the entanglement 
sudden death (ESD) phenomenon In most cases, in 
order for a quantum information task to be successful, 
the ESD time scale ought to be much longer than the 
computation or operation time scale. Therefore an open 
challenge is to determine and control the timing of ESD. 

A number of successful schemes, such as quantum error 
correction [f| , feedback control [1] , dynamical decoupling 
0, etc, have been proposed during the last two decades 
in order to control and suppress the decay of local co- 
herence. However, the decay of non-local coherence, i.e., 
entanglement, can be very different and inconsistent with 
local decoherence 0. Consequently, several recent stud- 
ies 0-11 1 based on different schemes have revealed that 
while these proposals may delay the occurrence of ESD 
in some cases, they may also counter-intuitively acceler 



ate or even induce ESD. Decoherence free subspace [12 1 
is another category of strategy. However, it would ei- 
ther cost an increased number of non-computing qubits 
to create certain symmetry, or some serious engineering 
of the system or reservoir with external contacts (equiv- 
alent to dynamical decoupling [l3|) that might bring in 
additional source of noise to neighboring qubits. A more 
reliable, simpler, and systematic method to delay or even 
avoid ESD would be desirable. 

Most of the prior treatments are managed during the 
decoherence and/or disentanglement process. An inter- 
esting question is whether ESD can be pre-managed be- 
fore it is taking place, i.e., during the state preparation 
process. Recently we have given a trial answer to this 
question for a type of simple form two-qubit states. We 
showed that ESD properties are determined by a couple 



of initial parameters, i.e., ent ang lement and double ex- 
citation probability or purity [l4j . It is also shown that 
ESD and initial pure states of a two-qubit system^ under 
identical noisy channels) are statistically related [l5j • For 
practical situations, the initial state may not be simple 
and pure, and the reservoirs may not be identical. Then 
a key question arises: is it still possible to manage ESD 
with only a few controllable initial parameters. 

In this paper we extend the analysis to arbitrary two- 
qubit mixed states, and treat decoherence via local differ- 
ent amplitude-damping dissipations. We decompose the 
mixed state into its "X" and "O" components ljl 17 1, 
and first analyze the X-form component. This is a legal 
density matrix by itself and its entanglement properties 
were first studied by Yu and Eberly [l8j. We show that 
the exact ESD-free, ESD-tolerable and ESD-no-go phases 
of the X-form matrix are explicitly determined, and thus 
manageable, through several of its initial parameters, i.e., 
entanglement and excitation probabilities. 

We then apply the static entanglement boundary re- 
lation for X-form matrices, developed recently by Ma et 
al. [IH and by Hashemi Rafsanjani and Agarwal fl7j j . 
to the dynamical evolution of the state. They show that 
the ESD phase critical conditions of the X-form compo- 
nent simply provide lower bounds for those of the entire 
state X+O. Therefore the control of different ESD prop- 
erties of a complicated state can be realized sufficiently 
through just a few initial parameters of its X-form com- 
ponent. Our results may provide guidance to experimen- 
tal preparations of robust entangled states. 



INITIAL STATE AND TIME EVOLUTION 

In practice, state preparations will inevitably leave the 
system in mixed states. Thus we consider the most gen- 
eral form of a two-qubit (A and B) initial state 
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which is expressed in the standard basis: |e)A|e)s, 
\ e )A\g)B, \g)A\e)B, \9)a\9)b- Here |e) and \g) are usually 
interpreted as excited and ground states of qubits A, B. 
It is obvious that such a matrix can always be decom- 
posed into a sum of X-form and O-form matrices, i.e., 
p(0) = p x (0) + O(0). Here 
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has the X shape matrix elements [18| . and O(0) contains 
the remaining matrix elements that form an O shape. 
According to Refs. 16|,|l7|, the entanglement in terms of 
concurrence [19] for the entire state is no less than that 
of its X-form component, i.e., C[p(0)} > C[px(0)]. 

After the state preparation, the two qubits will in- 
evitably interact with the environment and begin the 
process of dissipation. Here we consider the case when 
qubits A, B are subject independently to two different 
amplitude-damping channels, which may describe many 
physical realizations, e.g., spontaneous emission, spin 
chain Heisenberg interactions [2(j, etc. The evolution 
of the initial state can be described by Kraus operator 
representations [ljj, i.e., 



(3) 



where the four joint Kraus operators Ki(t), i = 1, 2, 3, 4, 
satisfying the relation Yli=i Ki(t)Kj(t) = 1, correspond 
to four possible tensor products of individual Kraus oper- 
ators of qubits A, B : i.e., Ki(t) = ® K^, p, v = 0, 1. 
For the amplitude-damping mechanism one has 
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(4) 



where p a is the excitation-transfer probability that grows 
from to 1 irreversibly, and q a = 1 — p a denotes the 
probability that qubit A remains in its excited state. The 
Kraus operators for qubit B are characterized similarly 
as (U]) with decay probability q b . 



ESD PHASES FOR X-MATRIX 

We first consider the dynamics of the X-form compo- 
nent ([2]), which can be regarded as an initial state by it- 
self, and analyze its ESD properties. Later we will extend 
the static relation C[p(0)] > C[px(0)] to the dynamical 
evolution of the entire arbitrary state ([lj, and show its 
ESD properties are manageable through the controllable 
parameters of the X-form component. 

The degree of entanglement, i.e., concurrence, between 
the two qubits for the X-form initial state @ takes a very 




FIG. 1: Schematic illustration of two types of two-qubit ma- 
trices. For convenience the height of the cuboid represent the 
modulus of each matrix element. Plots (a) and (b) illustrate 
the <3> type (Q<s> > Q<s) and type (Q* > Q$) of two-qubit 
matrices respectively. 



compact form [16 



C — max{Q$,(5*,0}, 



(5) 



where 



Q<S> = 2(|/5 14 | - v/p22/>33), Qf = 2(1/923 | - \fpHpM)- (6) 

Obviously, C = when Q$ < and < 0. This is 
not an interesting case, so we will consider the cases when 
one of the Qs is positive (matrix positivity does not allow 
both Qs to be positive simultaneously). Then the initial 
concurrence is expressed conveniently as Co = for 
Q$ > Q*, or Co = Q* when > Q$. Accordingly, we 
define the X-form matrix as well as the entire matrix ([T]) 
as a $ matrix if Q$ > Q^, oraf matrix otherwise. An 
imperfect preparation of a $-like Bell state is a practical 
example of generating such a $ matrix. Fig. [T] gives a 
schematic illustration of both types of matrices. 

By applying the Kraus operators ([!]) to the initial state 
@, one will obtain the time dependent state as 



Px{t) = 



{ pu(t) pu(t)\ 

P22{t) P2 3 (t) 

P 32(t) P33 (t) 

\ P 4l(t) pn{t) J 



(7) 



where pn(t) = pnq a qb, = piiq a Pb+P22q a , P33W = 

PllPaqb + P33Qb, Pu(t) = PllPaPb + P22Pa + P33Pb + PiA, 

Pu(t) = p|i(t) = puy/q^qb and p 23 (t) = p3 2 (t) = 
P23y/q a qb- The state retains an X form [18| . and con- 
sequently the time dependent concurrence can be simply 
obtained as C(t) = max{<5$(t), Q^(t), 0}, with 



Q<s,(t) = 2y/q a qb{\pu\ 
= 2^q a q b {\p23\ 



P22P33), (8) 
P11P44), (9) 



and 



A(t) = pll(pllPaPb + P22Pa + P33Pb) ■ 



(10) 



The irreversible damping process indicates that the prod- 
uct q a qb is decreasing and A(t) is increasing with time. 
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This guarantees the decay of both Q$(t) and Q*(t). 
Thus one has the time dependent concurrence, C(t) = 
max{<2$(i),0} or C(t) = max{Q$(t), 0} for $ or * type 
of initial matrices respectively. 



$ type of X-matrix 

Now we turn to the central task of this paper, and 
analyze the conditions of state ^ to avoid ESD, as well 
as the ESD onset conditions if it were to occur. We set 
Xij = y/PuPjj and yij = p u + p r] with i,j = 1, 2, 3,4 for 
the convenience of our following notations. 

For the 4> matrix situation, i.e., C'o = Q$ > , the 
matrix positivity requires Q$> /2 < Xu — X23- This can be 
transformed in terms of the double excitation probability 
D = p n , and leads to the relation: D% in < D < D£ ax , 
where the physical boundaries are given as 



2£>S 



1 - 2/23 - V (1 - 2/23) ~ (<9* + 2x 23 ) , 



2Dr x = 1-2/23 + V( 1 -y23) 2 -(Q* + 2a; 2 3) 2 (ll) 

Fig. [2] illustrates the behaviors of L>g lax and L>| lin with 
respect to Q$> by the upper red and lower blue solid lines 
respectively. The two boundary lines connect when D = 
(1 — 2/23 )/2 and Q$ takes its maximum value Q$ ax = 
1 — 2/23 — 2x23 as illustrated by point 1 in plot (a). 
From the quantity Q$(/) given in ([8]). the condition to 



have ESD is whenever A(£) 



! "23 



> I P14 1 . There is no 



ESD at any time if A max (t)+x% 3 < |pi4| 2 , where A max (Z) 
is the maximum value of A(t) and is achieved when p a , 
Pb take their maximum value 1 (usually corresponds to 
t = 00). This ESD-frec condition can be re-expressed in 
terms D bounded by its critical value D^, i.e., 



D < D 



viz 



Q% +4z 23 Q$ -2/23 /2. (12) 



Interestingly, D is only controlled by the system initial 
parameters P22, P33- So here ESD cannot be avoided 
by adjusting the strength of reservoir interactions, e.g., 
the decay rate of q a and qb with respect to time. 

Fig. [5] (a) is an ESD susceptibility phase diagram of 
the $ type of X-form initial state @ with respect to Q$ 
(Co) and double excitation D with the single excitations 
P22, P33 treated as fixed parameters. The black solid 
line is the ESD-free critical value and straightfor- 
ward calculations show that it can only have cross points 
with the minimum physical boundary -D™ 11 when Q$ = 
A - 2x 23 , and 2D = y/yf 3 - 2x\ z + A - y 23 , 

and 



\f 7 2x 



23 



23' 



Here 



with 2A = 1 - 27/23 ± y/(l ~ 22/ 23 ) 2 - 8a 
+ correspond to points 2 and 3 respectively as illustrated 
in the plot. From the condition (fl2|) . one identifies the 
region below in Fig. [5] (a) is the ESD-free phase and 




FIG. 2: ESD phase diagrams with respect to the double ex- 
citation D and X-matrix entanglement Qc> (Co) for $ type 
of initial states. Plot (a) illustrates the ESD-susceptible 
(shaded) and ESD-free (unshaded) phases separated by the 
black solid critical line for X-form initial state ((2}. Plot 
(b) is plotted on top of (a), and it illustrates the ESD time 
phases in the perspective of traffic-light signals. Dashed lines 
(i), (ii), (iii) represent the ESD-tolerable critical value D% for 
r = 7a/7b = 1,3, 10 respectively. The green zone below I)£ 
is the ESD-free (safety) phase. For r — 1, the ESD-tolerable 
(caution) phase is the yellow zone between and dashed line 
(i), and the ESD-no-go (danger) phase is the red zone above 
line (i). The vertical line illustrates Co = Q* = 0.8. For 
the generic state {TJ, t ne corresponding ESD-free and ESD- 
tolerable critical lines will move up a little as schematically 
shown by lower and upper dotted lines respectively. 



everywhere above is the ESD-susceptible phase. Obvi- 
ously, the double excitation needs to satisfy the relation 
£,min <D<D{ to guarantee an ESD free zone. 

We further analyze the ESD onset relation Q<s>(t) = 0. 
One can consider a practical situation that quantum com- 
putation tasks can tolerate ESD if it happens only after a 
certain time point r, which is much longer than the com- 
putation and operation time scales, and any ESD that 
happens earlier than that may be harmful and is consid- 
ered as a no-go region. Then the ESD-tolerable condition 
is given as t E sD > r or vu^esd) > Pfe(r), where k = a,b 
and tESD denotes ESD onset time. This relation can be 
expressed in terms of the double excitation D bounded 
by its critical value i.e., 



D< D% = 



P33 

Pair) 

P33 

Pb{r) Pair)' 



I P22 

\Pbir) 
P22 



(gg + 2x 23 ) ; 
Pa(r)p fc (r) 



(13) 



One can show dD^/dpkir) < 0, i.e., is a decreasing 
function of the excitation-transfer probabilities Pkir). 

For demonstration, one can reasonably assume the 
channel dissipation as an exponential decay process, i.e., 
Pk(t) = 1 — e _7fet , where jk (fc = o,,b) are the decay 
constants for qubits A, B respectively. Then it follows 
straightforwardly that D\ is also a decreasing function 
of 7^. Fig. H](b) illustrates the behavior of when one 
fixes channel 7 Q and varies the other channel 7&. The 
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dashed lines (i), (ii), (iii) represent the cases when the 
ratio r = j a /"fb = 1,3, 10 respectively. Here r is taken 
to be the decay lifetime l/j a - Note that 7 a ,7b = are 
singularities of D$ and they represent the relatively less 
practical single-channel dissipations. Their ESD-free and 
ESD-tolerable critical conditions can be obtained easily 
from (|SJ) and will not be addressed here. 

The ESD-susceptible phase in Fig. [5] (a) is then further 
separated into two phases. In Fig. [2] (b) the green zone 
is exactly the ESD-free phase in plot (a), and is consid- 
ered as the safety phase where entanglement decays to 
zero only at tESD = oo. For r = 1, the ESD-tolerable 
(caution) phase is the yellow zone where ESD occurs only 
after the tolerable time scale, i.e., tssD > t, and the red 
zone indicates a fast ESD that happens before r, which 
is the ESD-no-go (danger) phase that state preparations 
should avoid. Obviously, this ESD-tolerable region will 
increase with the increase of the ratio r. 



$ type of X-matrix 

For the *Sf matrix situation, i.e., Co = > Q<s>. The 
restriction of can be achieved as < D < D™ x , 
where the maximum physical boundary 

(14) 

is illustrated in Fig. ^ by the red solid line as a function 
of Qty. Under this boundary (j 14[) . takes its maximum 
value Q£ ax = 2^/p 22 (l-y 24 ) when D = 0. 

From (j9)) the condition to have ESD is whenever A(£) + 
x \i > |yC>23 1 2 - Then the ESD-free condition is achieved 
as A max (t) + x\ A < 1/923 1 2 , which again doesn't depend 
on the strength of reservoir interactions and can be re- 
expressed in terms of the double excitation probability 
bounded by its critical value D^,, i.e., 

D<d{ = Q%/A(1-^) 2 . (15) 

Fig- 131 ( a ) is an ESD susceptibility phase diagram of 
the "J type of X-form initial state (HJ with respect to 
Q* (Co) and double excitation D with the zero excita- 
tion P44 treated as a fixed parameter. The black solid 
line is the ESD-free critical value Dq, that separates the 
ESD-free phase from the ESD-susceptible phase, and it 
crosses the maximum boundary line D^ x at a single 
point (illustrated by point 1 in the plot) when =2(1— 

Vm)VP22 (y24 - yfJ/Vs and d = p 22 (y 24 - yl 4 )/ s , 

with S = (,022 + \fPil) 2 + P22(l - \fp~ii) 2 - Interestingly, 
Eq. (p~5)) and Fig. [3] (a) show that when the double exci- 
tation vanishes, i.e., D = 0, the *f? type of initial state is 
ESD free for any amount of (Co). 

From the ESD onset condition Q^(t) — 0, one can 
again re-express the ESD-tolerable condition tESD > t 




FIG. 3: ESD phase diagrams with respect to the initial dou- 
ble excitation D and X-matrix entanglement Q* (Co) for ^ 
type of initial states. Plot (a) illustrates the ESD-susceptible 
(shaded) and ESD-free (unshaded) phases separated by the 
black solid critical line for X-form initial state ([2]). Plot 
(b) is a re-plot of (a), and it illustrates the ESD onset time 
phases in the perspective of traffic-light signals. Dashed lines 
(i), (ii), (iii) represent the ESD-tolerable critical value D% for 
r — 7a/7b = 1,3, 10 respectively. The green zone below the 
black solid line is the ESD-free (safety) phase. For r = 1, 
the ESD-tolerable (caution) phase is the yellow zone between 
D-f, and dashed line (i), and the ESD-no-go (danger) phase is 
the red zone above (i). For the generic state (fTJ, the corre- 
sponding ESD-free and ESD-tolerable critical lines will move 
up a little as schematically shown by lower and upper dotted 
lines respectively. 

in terms of the double excitation probability as D < D\, . 
Here D\, is the critical value that can be obtained ana- 
lytically by solving for D from the equation 

Q|/4 = Pb {T)[p a {T)-l]D 2 -Q^^J5 

+ [P22Po.(t) +Pb(r) - U2iPb(r)]D. (16) 

The analytic dependence of D\, on the parameters 
Q*, P22, P44, Pa(r), and Pbij) is a very complicated 
expression. With the assumptions, Pk(t) = 1 — e~ 7fct 
and t = l/7a, the explicit behavior of D\, is illustrated 
in Fig. [3] (b) by the dashed lines (i), (ii), and (iii) for 
r = la Mb = 1,3,10 respectively. In the plot the green 
zone is exactly the ESD-free zone in plot (a), and it is 
the safety phase with ESD happens at time infinity. For 
r = 1, the ESD-tolerable (caution) phase is the yellow 
zone with ESD happens after the tolerable time r, and 
the ESD-no-go (danger) phase is the red zone with ESD 
happens before r. Again, the ESD-tolerable region in- 
creases with the increase of the ratio r. 



ESD PHASES FOR ARBITRARY MATRICES 

Now we return to the generic initial state ([1]), and note 
that the evolutions of px(0) and O(0) are independent 
of each other. That is, p x (t) = T,t=i K i{t)Px(0)Kj{t) 
remains in an X form as shown in l[7]). and 0(t) — 
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T,i=i K i(t)O(0)K^(t) retains the O form, i.e., 



O(i) 



( Pl2 (t) p 13 (t) \ 

P2l(t) p 24 (t) 

Psi(t) p 34 (t) 

/>«(*) P43(i) 



(17) 



where p u (t) = p* 21 (t) = pi 2 q ay /qb, = pt\{t) = 

PuQby^, P2i(t) = P* 42 {t) = P24y / ql + Pi3Pb^/qI, and 

P34(t) = p* 43 (t) = p 3 4y/^b + P\lVarM>- 

Then the analysis in Refs. [16J, LLZ| for static matri- 



ces can be extended to the time dependent state p(t) = 
Px(t) + 0(t). One immediately notes that the time de- 
pendent concurrence of p{t) is always no less than that 
of its X-form component, i.e., C[p(t)] > C[px(t)} at any 
time t. Consequently, the arbitrary state p(0) will always 
reach ESD no earlier than px(0). This is to say that the 
exact ESD- free critical conditions (fT"2|) and (fT5j) for px (0) 
will simply provide lower bounds for those of the entire 
state p(0). Schematic illustrations of those ESD-free crit- 
ical lines for p(0) are provided in Fig. [2] (b) and Fig. [3] 
(b) by the lower dotted lines, which are guaranteed to be 
no lower than the corresponding exact ESD-free critical 
lines and for px{0)- The same argument is also 
valid for the ESD-tolerable critical conditions, which are 
schematically illustrated by the upper dotted lines in the 
corresponding plots for r — Ja/jb = 1- Then the ESD- 
free, ESD-tolerable and ESD-no-go phases separated by 
the two dotted critical lines can be identified schemati- 
cally for the generic initial state (jTJ) . 

The key information here is that those exact critical 
conditions for the X-form matrix ^) arc sufficient con- 
ditions to guarantee the arbitrary state ([1| to be ESD 
free and ESD tolerable respectively. In this sense, the 
pre- management of ESD properties for the complicated 
general state p(0) can be simply realized by managing 
just a few parameters that explicitly determine the ESD 
properties of its X-form component px(0)- 

Experimental preparations of initial parameters may 
not always be precise, e.g., to produce a desired amount 
of initial entanglement Co is a challenge. However, our 
ESD phase analysis shows that one is still able to control 
the ESD properties within a range of these parameters. 
To be more specific, let us take the $ type of X-form com- 
ponent px(fi) as an example. Suppose the experimenter 
is able to prepare initial entanglement within the range 
0.80 < C*n < 0.96. 



tations, e.g., p 22 
shows that £>™ ax 
and D™ in 



This requires very small single exci- 
P33 — 0.01. Some simple arithmetic 
0.76, D% ~ 0.63 (r=l), L>£ ~ 0.40, 
0.22 for the lower bound case Co = 0.8, illus- 
trated by the vertical dashed line in Fig. [2] (b) . Then the 
ESD-free, ESD-tolerable phases can be realized when the 
experimenter controls the double excitation in the ranges 
0.22 < D < 0.40, 0.40 < D < 0.63 respectively. 

The monotonic increasing nature (with respect to Co) 
of critical values and D*^ indicates that the preferred 



ranges of D for Co = 0.8 will guarantee the corresponding 
ESD properties for states with Co > 0.8. Moreover, those 
ranges are also guaranteed to assure the corresponding 
ESD properties for the generic state p(0). 



CONCLUSION 

In summary we have studied possible pre-management 
of non-local decoherence, i.e., ESD susceptibility and 
timing, with limited resources of initial conditions. This 
is an open question not only from the view of funda- 
mental entanglement theory, but also from the practical 
aspect of preparing robust entangled states. 

It is easy to imagine that ESD properties may relate, 
at least partially, to initial conditions. However, to pre- 
manage them for arbitrarily complicated states with just 
a few initial parameters is not so intuitive. Our results 
show that it is possible to realize such a pre-management 
through initial entanglement and excitation probabilities 
for a generic two-qubit state that undergoes local differ- 
ent amplitude damping dissipations. These parameters 
just amount to a set of management tools that control 
the critical conditions for ESD-free, ESD-tolerable, and 
ESD-no-go phases of the state. In addition, these phases 
allow a well defined amount of margin for errors in prac- 
tical imprecise state preparations. 

We remark that the cases when the environmental 
noises are of phase-damping or mixed-damping mecha- 
nisms will show similar ESD phases. These will simply 
enlarge the picture of pre-manageable ESD phenomenon 
and will be addressed elsewhere. 
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